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Abstract 

Computational and experimental results on the thermally-induced magnetization reversal in 
single-domain magnetic nanoparticles are reported. The simulations are based on the direct inte- 
gration of the Fokker-Planck equation that governs the dynamics of the magnetic moment associ- 
ated with the nanoparticles. A mean field approximation is used to account for the influence of the 
dipolar interaction between nanoparticles. It is shown that the interactions can either speed up or 
slow down the reversal process, depending on the angle between the external magnetic field and 
the axis of easy magnetization. The numerical results are in good agreement with experimental 
measurements on cobalt-platinum nanoparticles. 

PACS numbers: 75.20.-g; 75.60.Jk; 75.75.Jn 
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I. INTRODUCTION 



Single-domain magnetic nanoparticles constitute an attractive system for fundamental 
research as well as for advanced technological applications. The use of single-domain mag- 
netic nanoparticles is expected to increase the data storage density to several petabit/inch 2 
(10 15 cm -2 ) [l, 2] in the near future. However, when the size of the nanoparticles is reduced, 
the superparamagnetic regime can be attained and the magnetization fluctuates under the 
effect of thermal excitations Jsj]. This effect is a major drawback for technological devel- 
opments and it is therefore important to investigate and understand the thermally-induced 
magnetization reversal in these systems. 

In isolated single-domain magnetic nanoparticles, the magnetization reversal by thermal 
activation is well described by the Neel-Brown model % |5[. According to this model, the 
thermal fluctuations cause the magnetic moment to undergo a Brownian-like motion about 
the axis of easy magnetization, with a finite probability to flip from one equilibrium direction 
to another. From an energetic point of view, the two minima associated to the equilibrium 
positions are separated by a barrier due to the magneto-crystalline and shape anisotropics. 
The corresponding Arrhenius-type superparamagnetic relaxation time can be written as 
t = TQexp(AEfkBT), where AE is the energy barrier between the two easy directions of 
the magnetization, ks is the Boltzmann constant, and T is the temperature. The typical 
time To is not well known experimentally and is estimated to be of the order 10~ 10 — 10~ 12 s 
for magnetic nanoparticles. 

The Arrhenius exponential law describes well the behavior of isolated (i.e., non- 
interacting) nanoparticles. The effect of dipole-dipole interactions on the relaxation time 
and, more generally, on the reversal process has been studied in several works, both theoret- 
ical js, 7| and experimental js-ll]. In spite of significant progress, it is still a controversial 
issue, as opposite dynamical switching behaviours have been reported in the literature jfj]. 
This is of course a major issue for the development of smaller and faster switching memory 
devices. 

In this work, we present numerical calculations of the relaxation times for the magneti- 
zation reversal in isolated and interacting single-domain ferromagnetic nanoparticles with 
uniaxial anisotropy. Further, in order to validate our theoretical approach, we have studied 
experimentally the dynamics of close packed cobalt-platinum core-shell nanoparticles which 
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can be considered as a model system owing to their spherical shape and small size dispersion 
(less than 5%). The relaxation time is determined by performing magneto-optical measure- 
ments using femtosecond laser pulses. As shown in the following sections, our theoretical 
model is in good agreement with the experimental results. 



II. MODEL 



Micromagnetic simulations of thermally activated magnetic systems are generally per- 
formed using either the Langevin Dynamics (LD) [l2], which is based on the direct integra- 
tion of the stochastic Landau-Lifshitz-Gilbert equation j^J, or the so-called 'time- quantified' 
Monte Carlo method (TQMC) [13], which is a generalization of standard Monte Carlo tech- 
niques to the time-dependent regime. The TQMC method was recently extended to the 
case of correlated particles with nearest-neighbors exchange interactions [14|, although this 
approach was limited to small systems (10 x 10 square lattice) with short-range interactions. 

Adapting the above methods to the case of long-range dipolar interactions would re- 
quire prohibitive computational times and memory storage, particularly when the number 
of nanoparticles is large. In order to circumvent this problem, we have developed an approach 
based on the direct integration of the Brown- Fokker- Planck equation {5J], which describes 
the time evolution of the probability distribution W(9,(f>,t) of the magnetic moment of a 
nanoparticle. 9 and <ft are the polar angles determining the orientation of the magnetization 
vector m with respect to the axis of quantification z (0 < 9 < 7r) and its projection on the 
equatorial plane (0 < <fi < 2tt). In the following, we consider the case of axial symmetry, 
for which all quantities do not depend on 0. Dipole-dipole interactions are treated in the 
framework of the mean field approximation. 

For an isolated nanoparticle in an external magnetic field H aligned along the axis of 
easy magnetization, the Brown- Fokker-Planck equation reads as jsj: 

dW 1 d 



dt sin 9 09 



. n f,,dE Txr ,,dW\ 



(1) 



k B T (l + a^VMs' y> 
where V is the volume of the particle, a is the Gilbert damping constant, 
7o=1.76 xlO 11 T _1 s _1 is the gyromagnetic factor, and M s is the magnetization at satu- 
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ration. The total energy of a nanoparticle is given by E{6) = KV sin 2 6 — /io m ■ Ho, where 
K is the anisotropy constant and /io is the vacuum magnetic permeability. 

It is useful to rewrite Eq. ([T]) using dimensionless quantities, by introducing the normal- 
ized time variable t = t ^"^f^- and defining x = cos 6, with —1 < x < 1. The resulting 
equation can be put in the form of an advection-diffusion equation: 

dW 

where D(x) = 1—x 2 , U(x) = (a; 2 — l)(Ax + B), and A and B are two dimensionless constants 



representing respectively the anisotropy energy and the energy due to the external field: 

A= 2KV B = t ^MJU 

k B T k B T y ' 




FIG. 1: (color online). A sketch of the initial (red line) and final (blue line) distributions of the 
moments with respect to the energy barrier (green line). 

We solved Eq. ([3]) using a finite difference technique with boundary conditions chosen 
so that the total probability J. W(x,t)dx remains constant in time. Initially, the magnetic 
moments are set near the local equilibrium with higher energy, which is situated at x — — 1 
{9 = 7r), so that the probability distribution W{x, 0) is given by a half-Gaussian with 
maximum at x — — 1 (red curve on Fig. [T]). Physically, this corresponds to a metaequilibrium 
with a magnetic moment oriented in the direction opposite to the external magnetic field. 
Our procedure consists in computing the average time necessary for the magnetic moment 
to reverse its orientation and align with the external field. This configuration minimizes the 
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energy for all angles 6 and therefore constitutes an accurate description of the equilibrium, 
given by the Boltzmann law Wfi na i = const, x e ~ E l k BT ( D m e curve on Fig. [1]). 

In order to obtain the relaxation time, we compute the integral I-(t) = j^_ l W{x,t)dx, 
which represents the probability of finding the magnetic moment on the lower part of the 
equatorial plane, and then we fit this quantity with a decaying exponential. This procedure 
is repeated for different values of the temperature. For comparison with experimental re- 
sults, we will apply our numerical simulations to the case of cobalt nanoparticles, which are 
ferromagnetic at low temperature (typically, the superparamagnetic blocking temperature 
is ~ 100 K for nanoparticles with a diameter less than 10 nm). 
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FIG. 2: (color online). Relaxation times for isolated nanoparticles (circles) and interacting nanopar- 
ticles for different interparticle distances: d = 9.6 nm (squares), d = 11.8 nm (triangles), and 
d = 16.6 nm (stars). The occupation probability is p = 0.5. The damping constant a = 1. 




— * — interacting: d=16.6 nm 
— A — interacting: d=l 1.8 nm 
— ■ — interacting: d= 9.6 nm 



III. NUMERICAL RESULTS 

First, we focus on the case of non- interacting (isolated) particles. We considered spherical 
Co nanoparticles with a diameter of 9.5 nm, an anisotropy constant K = 4.2 x 10 5 Jm -3 , 
and a saturation magnetic moment m s = V r M s =6.45xlO -19 JT _1 . The external magnetic 
field is H Q = 1.59 x 10 5 Am- 1 . 

We computed the relaxation time r as a function of the density barrier AE = KV(1 — 
fi H M s /2K) 2 . The results for isolated nanoparticles are plotted in Fig. [2] (black circles) 
and are in agreement with Brown's expression (dashed line) at low temperatures, with a 10% 
accuracy. At higher temperatures, the numerical solution departs from Brown's expression 
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and the computed relaxation time tends to zero. These results are consistent with those 
obtained using different numerical techniques such as time-quantified Monte Carlo [13] and 
thus fully validate our approach. 

Let us now focus on the more challenging problem of magnetic nanoparticles interacting 
via the magnetic dipolar interaction. The case of two particles was considered by Rode 



et al. 



15]. The generalization to many-particle systems is far from trivial jo, |7J, as the 



computational complexity increases as the square of the number of particles N 2 . This can 
be reduced to N\nN by making use of fast Fourier transforms 16, but still requires 
considerable computing time for large systems. 



Here, we use an alternative approach based on the mean field approximation [18l-l20|. 
whereby the effect of the dipolar interaction on one particle exerted by all the others is 
expressed as a self-consistent magnetic field. This approximation neglects higher-order fluc- 
tuating correlations between the magnetic moments, which are a further source of random 
Brownian motion beyond the usual fluctuations induced by thermal effects. In the regimes 
considered here, which operate at relatively high temperatures, the thermal fluctuations are 
likely to be dominant so that the mean field approximation should be accurate enough. 

We consider a system of interacting nanoparticles distributed over a spatial lattice. At 
sufficiently low temperature (or high energy barrier), the magnetic moments fluctuate in 
the vicinity of the positive and negative directions of the z axis. In this case, the self- 
consistent field is also aligned along the z direction and can be written as = Hu(t)e z = 
8pSd~ 3 m z (t)e z , where S is a constant that depends on the geometry of the lattice, < p < 1 
is the probability of occupation of the sites, d is the center-to-center interparticle distance, 
and m z (t) is the average z component of the total magnetic moment of the system: 

/'7T 

m z (t) = m s cos 6W (6, t)d(cos 6). (5) 
Jo 



For a two-dimensional lattice [18(, the mean field model yields a value S ~ —1.129, where 
the negative sign signals that the dipolar interaction is antiferromagnetic. 

In the mean field approximation, the nanoparticles are considered to be independent, 
with all particles experiencing the effective dipolar mean field Hd(£) as an external magnetic 
field. Thus, the total energy becomes E{6) = i^\^sin 2 6 — fiom ■ (Ho + Ho). This entails a 
modification of U(x) in Eq. (J3]), which becomes U(x,t) = (x 2 — l)[Ax + B + C(t)}, where 
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C(t) is a term describing the strength of the dipolar interaction: 

We emphasize that, within this approach, the relevant Fokker-Planck equation ([3]) becomes 
nonlinear, because the dipolar term depends on W through Eq. (jSJ). 

The numerical results for interacting particles are shown in Fig. 12} for the case with 
p = 0.5 and three values of the interparticle distance d. The Gilbert damping constant 
is taken to be a = 1, a value that is consistent with recent measurements on small Co 



nanoparticles (21|. It is found that the relaxation time decreases with decreasing interparticle 
distance, in agreement with previous experimental measurements [8j]. Beyond a certain 
distance (d w 16 nm in the case of Fig. |2j), the dipolar interaction becomes negligible and 
the result for the corresponding isolated nanoparticles is retrieved. In contrast, the relaxation 
time decreases with increasing site occupation probability p. In summary, the magnetization 
reversal process is accelerated for short interparticle distances and large concentrations. 

In all the preceding results, the external magnetic field H was set parallel to the axis 
of easy magnetization. When the field makes an angle \& with respect to the easy axis, the 
energy profile becomes: 

E{6) = KV sin 2 9 - ^ m s H Q cos(# - *) - /i m • H D . (7) 



Coffey et al. [22[ found that the relaxation time in isolated nanoparticles varies with the 
angle and has a minimum at \& = tt/A and a maximum at \& = tt/2. Our mean field 

n 

model [19[ assumes that the self-consistent dipolar field is always aligned along the axis of 
easy magnetization, which is a reasonable assumption only when the external field is either 
parallel = 0) or antiparallel = n) to z. 
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FIG. 3: (color online). Relaxation time r vs the normalized inverse temperature, for an external 
field antiparallel to the direction of easy magnetization (vf = 180°). The inset shows the same 
result for ^ = 0°. 

We have studied the relaxation time for nanoparticles with the following physical param- 
eters: K = 4.2 x 10 5 Jm" 3 , m s = 8 x KT^JT -1 , V = 5.6 x lfr 26 m 3 , d = 8.6 x l(T 9 m, 
S = —1.129 and p = 1. The results are shown in Fig. [3] for different values of the lattice 
temperature. When the external field is parallel to the z axis = 0), the effect of the 
dipolar interactions is to shorten the relaxation time, i.e., to accelerate the reversal of the 
magnetization. In contrast, for the antiparallel configuration = ir), the relaxation time 
is longer when the dipolar interactions are taken into account, at least for temperatures 
below a certain threshold (here the threshold is situated near KV/ksT = 3, corresponding 
to a temperature T = 567 K). This effect can be understood by noticing that the dipolar 
field either adds or subtracts to the external magnetic field, depending on the orientation 
of the latter. Nevertheless, since the dipolar field - and thus the energy barrier - depends 
self-consistently on W(6,t), nonlinear effects might play an important role in the reversal 
dynamics. Further work will be necessary to clarify these issues. 



IV. COMPARISON TO EXPERIMENTS 



In order to compare our theoretical results to experiments, we have measured the relax- 
ation times in core-shell cobalt-platinum nanoparticles. The CoPt core-shell nanoparticles 
are made by a redox transmetalation reaction. They have a spherical shape made of a cobalt 
core with an average diameter of 5 nm and a platinum shell with average thickness 1.5 nm 
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as measured by high resolution electron microscopy (JEOL 2100F microscope). The size 
dispersion is less than 5%. These nanoparticles are assembled into a compact bulk pellet 
by cold pressing under 160 Pa. Without any thermal process, the nanoparticles display a 
superparamagnetic behavior at room temperature with a blocking temperature of 66 ± 3 
K as measured by Zero Field Cooling/ Field Cooling in a SQUID apparatus. At 5 K, the 
nanoparticles are ferromagnetic with a coercive field H c = 7.56 x 10 3 Am -1 and a magneti- 
zation at saturation M s = 4 x 10 4 Am -1 . The inter-particle distance (Co-Co) is estimated 
to be 8.6 nm and the particle concentration is large enough so that in the simulations we 
can safely assume p — 1. 

To measure the relaxation times in core-shell cobalt-platinum nanoparticles we have per- 
formed time-resolved pump-probe magneto-optical Kerr measurements. The probe pulses 
are obtained from an amplified Ti: Sapphire laser functioning at a wavelength A = 800 nm 
and with a duration of 120 fs. The pump pulses with a wavelength A = 400 nm are obtained 
by frequency doubling in a (3 — BaB204 nonlinear crystal and have a duration of 150 fs. 
Upon heating the nanoparticles, the pump pulses also trigger a motion of precession of the 
magnetization that can directly be observed on the time resolved magneto-optical Kerr sig- 
nals. The experimental times r are obtained by fitting the relaxation of this laser-induced 
precession, and are determined for several values of an external static magnetic field (which 
is parallel to the direction of easy magnetization, i.e., ^ = 0). 

It must also be stressed that the laser pulse penetrates into the sample for a distance 
of the order of 10 — 15 nm, corresponding to just a few layers of nanoparticles. This pen- 
etration depth is much smaller than the sample dimensions, so that the excited region is 
approximately two-dimensional, in accordance with our theoretical assumptions. 

The important physical parameter that has to be determined is the lattice temperature. 
The nanoparticles are excited by a short laser pulse that heats up the electron gas, which 
subsequently cools down by exchanging energy with the lattice. After a few picoseconds, 
the electrons and the lattice are in thermal equilibrium at a temperature that is usually a 
few hundred degrees above the initial temperature. This equilibrium temperature can be 



estimated using a simple two-temperature model [23|, [24J, which for parameters relevant to 
the present experiment yields an equilibrium temperature around T = 520 K (inset of Fig. 

SD. 
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FIG. 4: (color online). Relaxation times for three values of the external magnetic field: experi- 
mental results (black squares); simulations including interactions (red circles); simulations without 
interactions (blue triangles). The simulations were carried out at a temperature T = 520 K. The 
inset shows the evolution of the electron (black solid line) and lattice (red dashed line) temperatures 
following laser excitation, obtained with a two-temperature model. 



This value of the temperature was used in our Fokker-Planck simulations of the magne- 
tization dynamics. The results are plotted in Fig. H] and show a good agreement between 
the experiments and the simulations that include the dipole-dipole interactions. It must be 
pointed out that the experimental relaxation times are subject to large errors due to the 
difficulty of determining precisely the Gilbert damping, so that an agreement between theory 
and experiment within a factor of two is already remarkable. Importantly, the simulations 
without interactions are clearly in much poorer agreement, for the absolute values as well as 
the general trend. 



V. CONCLUSION 



We developed a Fokker-Planck model that allows to simulate the magnetization dynamics 
in isolated and interacting single-domain magnetic nanoparticles. Isolated nanoparticles 
follow the Neel-Brown Arrhenius-like law at low temperatures (or, equivalently, large values 
of the energy barrier), but deviate from it at high enough temperatures. In the case of 
interacting nanoparticles, we made use of a mean-field approach, whereby the magnetic 
moment of each particle interacts with the mean dipolar field generated by all the others. 
Significant deviations from the Arrhenius law were observed, leading to a faster reversal 
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process in the presence of the magnetic dipolar interaction. The simulation results were 
in good agreement with experimental measurements performed on Co-Pt nanoparticles. 
Finally, we showed that, when the external field is antiparallel to the direction of easy 
magnetization, the effect of the dipolar field can be to slow down, rather than speed up, the 
reversal process. 
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